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Abstract: Systems of planar electrodes arranged on dielectric or piezoelectric layers are 
applied in numerous sensors and transducers. In this paper electrostatics of such electrode 
systems is presented and exploited in the analysis of distributed piezoelectric transducer 
dedicated to surface elastometry of biological tissues characterized by large Poisson 
modulus. The fundamental Matlab® code for analyzing complex planar multiperiodic 
electrode systems is also presented. 
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1. Introduction 

In many sensors, electric field is applied to investigate body by multiple electrodes which frequently 
can be considered periodic. This for example is the case of surface wave sensors of gas, actuators and 
linear ultrasonic motors utilizing planar metal strips as electrodes. Analysis of the field distribution and 
the electric property of electrode systems is usually necessary for the design and evaluation of the sensor 
parameters. The analysis of conducting strip is the general subject of this paper, and its application for 
elastometry is proposed as an example of its usefulness. 

Typically, elastic properties can be evaluated by measuring the ultrasonic wave time of flight over 
certain distance. In the case of tissue however, small dimension of tissue sample and considerable wave 
damping makes the task difficult. In this paper we show how to measure the wave velocity within one 
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ultrasonic transducer applying multiple strip electrodes on a piezoelectric layer applied to the tissue 
sample, in which case the detailed analysis of the strip system is indispensable because the frequency 
characteristic of the system, instead of the time of flight, is measured. 



2. A Template Electric Field 

The known identity [1] lays foundation for the presented analysis: 



r (l ~ E P " ( COS 9 ) COS (" + l) V 



'tt/2 sin' 4 9 



0<v <9, 



(cos v — cos 9)^ +l l 2 

0, 9 < V < 7T, 



(1) 



0 < 9 < it, Re{yii} < 1/2. The involved Legendre function P n (-) = P£{-) for \i — 0 has the following 
properties (for arbitrary real v and n integer) 



P—iz—l — Pui Pn{ %) — S n ( 1) P n (x), S n 



1, n > 0, 

1, 71 < 0, 



(2) 



Two equations can be casted from the above identity which can be easily interpreted in electrostatic 
terms (the second equation results from the first one after substitutions v — > v — n and $ — > i3 — n) 



/2)Kx 



]T P n (A)e- j ^ n+1 



S n P n (A)e-^ +1 ^ 



, \x\ < w, A = K w, 

V cos Kx — A 

0, w < \x\ < A/2, 



0. 



(3) 



\x\ < u>, 



-jS x V2 



V A — cos Kx 



, w < \x\ < A/2, 



where the variable d is replaced by Kx in order to accommodate the results to our purposes; 
A = cosA'w. 

Figure 1. Periodic system of strips with external cross-less connections within the strip cells 
including five strips; the arrangement used in the discussed sensor. 
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The electrostatic interpretation of the above complementary set of functions is the following. The first 
equation is the Fourier expansion of the normal electric induction on planar periodic perfectly conducing 
strips arranged along axis x with period A = 2n/K (K is the corresponding "wave-number" of strips), 
and width 2w. As known from electrostatics, the electric charge on strips equals the normal induction 
discontinuity across the strips (cf. Figure 1) that is square-root singular at the strip edges. The second 
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equation represent the Fourier expansion of the tangential electric field (similarly singular at the strip 
edges), naturally vanishing on the perfectly conducting strips (as opposite, the electric charge vanishes 
in the space between strips). These two expansions will be exploited as a template electric field in plane 
of planar system of periodic strips. Note that the Legendre functions involved in the left-hand side of the 
equations can be easily evaluated with help of efficient and robust numerical algorithm [2]. 

In order to make use of the above set of expansion, we need certain characterization of the electrostatic 
medium the strips to be embedded in. Consider the dielectric halfspace z > 0 of dielectric constant e. 
Applying the electric potential tp = exp(-jpx) on the plane z = 0 (which axis is considered normal 
to the plane of strips), resulting in tangential electric field E x = —V x tp = jpexp(-jpx), one easily 
evaluate that the excited normal electric induction (at z = 0+) satisfying the condition of vanishing field 
at z — y oo is (superscript '+' marks the considered field in the upper halfspace): 

Dt = -jeS p E+ (4) 

(neglecting the exponential terms), which is exactly the case of the corresponding harmonic components 
of the above two complementary Fourier series provided that the first equation is multiplied by e. The 
corresponding field components are P n (A) exp(— jnKx) and S n P n (A) exp(-jnKx), for D z and E x , 
respectively. For the other halfspace z < 0 with dielectric constant e and the field vanishing at z — y — oo, 
the corresponding equation to the above one is D~ = jeS p E~, which allows us to evaluate the electric 
charge distribution on the plane z = 0 between these two halfspaces as D = D+ — D~ = (1 + e/e)D+. 

3. Arbitrary Potentials on Strips 

The application of different electric potentials to subsequent strips breaks the electric field periodicity, 
hence the above Fourier expansion must be generalized into Bloch expansion of the planar electric field 
components by corresponding multiplication by a m exp(—j(r + mK)x) where r e (0, K) is constrained 
to one Brillouin zone for uniqueness reason, yielding (p„ = r + nK and series are simply rearranged; 
the superscript "+" is dropped, as well as subscript "x" at tangential field and z at normal induction): 

oo 

n=— oo m /c\ 
oo \P) 

E ( x ) = Y.3 S n-m a m P n . m (^)e~ jPnX 

n=—oo m 

(a rn are arbitrary constants and e e = e + e is the surface effective dielectric permittivity). This 
multiplication does not change the support domains of the normal induction D (which remains 
x E (—w, w) in the period A) nor x ^ (—w, w) for the tangential field E. 

Usually, the most interesting for applications are the strip potentials (which frequently are given) 
and the resulting strip charges or currents J = jujD, where to is the angular frequency of the 
applied harmonic potentials to strips Viexp(jtot) (I is the particular strip number). In the considered 
case characterized by Equation (5), both these quantities can be evaluated explicitly using Dougall's 
identity [1], which, applying Equation (2), can be transformed into two equations: 

= sin^vr y g^(A) = su^ 1 y (-1)"P B (A)_ 
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The / th strip potential V% is evaluated by simple integration of electric tangential field, and, knowing 
that this potential is constant over the entire strip, the integral can be evaluated at centers of strips, that 
is at x = I A (cf. Figure 1); it is for given value of r: 

Vj(r) = - J E(x)dx 

p—j(r+nK)x I 

= -3 £ a„g„-„P„-„(A) _ j{r - Jg^ = a ro K ( ""e-' rfA , (7) 
K (m| = lE^T% = -r^-irP-^rf-A). 

K n' + m + r/ K smwr/K 

n' 

Analogous integral of the normal induction D has constant values between strips (where D = 0), 
which allows us to evaluate this integral at the center of space between strips (that is at iA — A/2 or at 
iA + A/2 on both side of i th strip). The value of the i th strip charge is the difference of the integral values 
at these two points: at the space centers after and before the i th strip: 

Qi(r) = / D dx\ x={i+1/2 ) a - / D dx\ x=[i _ 1/2 )A = ea m Q^ e~ jriA , 

,A/2 ^ (8) 

/ J2 P n - m (A)e- j(r+nK)x dx = AP_ m . r/K (A). 

J -A/2 „ 



Q 



(m) 



-A/2 

Note that both Q and V are the spectral functions of r (moreover, a m may depend on r as well). They 
are, in fact, the Fourier transforms of discrete functions Qi and V\ of values taken at the strip centers 
x = /A. 



4. Strips on Layered Media 

If strips are placed on a layered substrate, the Equation (4) is no longer valid as the value of e depends 
on the component's wavenumber of the Bloch series, as it is shown in Appendix A for particular example 
of piezoelectric layer placed on a top of homogeneous elastic halfspace. The relation governing planar 
wavefield is: 

D(p) = -je(p)S p E(p), e(|p| > NK) w e e . (9) 

The fundamental feature of such system is that for large wavenumber value (say, for |r + nK\ > NK, in 
certain acceptable approximation), the surface effective permittivity reaches its constant limit e e , making 
Equation (4) valid for any n ^ [— 1 — A, N] (assuming r > 0 within the allowed limits), that is for 
infinite number of Bloch components in the expansion Equation (5). 

We apply the field expansion Equation (5) including summation over m with weight a m , which 
effectively means that each Bloch components of D and E is the sum like: 

oo 

D r (x) = D ne~ JPnX , D n = £f 5> m P n _ m (A), 

^=-°° m (10) 

E r {x) ^ E n € ^ , E n ^ ^ j$n— rnP^mPn— m(^) ) 

n=— oo m 

{p n = r + nK), where the appended subscript r marks the expansion dependence on spectral parameter 
r e (0, K). The summation limits over m, depending on N, is usually quite small and can be established 
by careful inspection of the equations discussed below. 
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These equations result from Equation (9) applied for any pair of harmonic amplitudes (D n , E n ) of the 
same harmonics exp(— jp n x) in the above expansion: 

- S n _ m e(p n )/e e ]P n _ m (A)a m = 0, (11) 

m 

for any n G (— oo, oo). Introducing e n = e(p n )/e e , the expression in brackets becomes: 1 — e n S n - m 
which is zero for |n| — > oo provided that m is finite. Noticing that e n = S n for n ^ [— 1 — N, N], 
according to Equation (9), it is easy to check that the discussed expression in bracket 1 — S n S n - m = 0 
for any n < —l — Norn > N provided that — 1 — N < m < N. This yields the condition for acceptable 
domain of m in Equation (10), which effectively yields the proper truncation of the discussed infinite 
system of Equations [3]. It is evident that choosing larger N and wider range of m, according the above 
rule, will yield only trivial solution to the additionally included unknowns a m . 

Note that the number of unknowns a m in Equation (10) is by one larger than the number of equations 
Equation (11), for \n\ < N. The required additional equation results from the Kirchhoff's laws applied 
to the system of strips; the simplest one is the condition on the strip potentials exploiting Equation (7). 

Evaluation of the strip voltages and currents (charges) requires integration of Equations (7, 8) over 
the spectral variable r, that is: 



K 

k&kh 

(12) 



Vl = Jc i ^ ®mV} m) e- jrlA dr = V k 

if 

t = jcuQ k = jue^ jT a m Ql m) e- irkA dr 



where to = 2%f is the angular frequency of the applied voltage to strips (/ is the frequency), 5^ is the 
Kronecker delta, V k are the given strip potentials, and l,k E (—oo, oo) are the strip numbers. It is evident 
from the first of the above equations that 

5> m (r)K (m) = ^e>' r,A . (13) 

m 

This is the equation which, appended to Equation (11), allows one to evaluate all unknowns a m which, 
substituted into the second equation presented above, finally yields [4] the admittance relation for simple 
periodic strips: 

I k =jue e Vi R(r)e-^ k - l)K dr/K, (14) 
Jo 

where I k is the current flowing to A;th strip of unitary length, resulting from the / th strip potential V\, and 

of \ Em«m P -m-r/if(A) . T 

R(r) = — — ; ; — — sm7r — . 

5. Simple Elastometric Sensor 

It is evident that the strip admittance depends on the piezoelectric and elastic property of the layered 
media on which the strips reside. In the case of biological tissue of large Poisson module (y < 1/2), 
the most important and difficult task [5] is the determination of the shear wave velocity v t of the tissue, 
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as the longitudinal wave velocity v\ can be easily measured. The discussion below presents a method of 
evaluation of the shear wave wavenumber k t = u/v t by means of strip admittance measured over certain 
frequency band. 

In typical cases, the spectral function R{r) is singular at the wavenumber r = fc R of Rayleigh 
wave (its value reduced to first Brillouin zone in the case of periodic system), but in order to make 
the measurements simpler in the considered case of elastometry, the layered system is chosen such that 
the Rayleigh-Lamb wave cannot propagate or it is not excited (see Appendix A for details). This feature 
allows one to perform simple numerical integration of Equation (14); below, the Fast Fourier transform 
is exploited for this task. 

Figure 2. Typical frequency dependence (a) of normal induction excited on piezoelectric 
layer residing on a tissue-like body, and (b) of corresponding signal S in infinite 
periodic system. 




The fundamental properties of the considered layered system are described by its effective dielectric 
permittivity e(r; iS) for chosen value of r and variable u, presented in Figure 2(a). This is the system 
response (in the complex amplitude of the normal induction) to the applied tangential electric wave-field 
on the system surface in the form exp j{ut — rx). One can easily notice that the effective generation of 
shear waves propagating along the system surface takes place at frequency f 0 = u Q / (2tt) of synchronism 
of the propagating wave with the delivered surface electric field spatial pattern: co 0 = rv t . It should be 
noted that similar phenomenon of generation of subsurface longitudinal waves takes place for much 
lower frequency rvi. The proposed and analyzed measurements concern, in fact, evaluation of e(r; uS) by 
means of the strip admittance dependence on angular frequency u. 

It is convenient to present first the results for infinite periodic system of 5-strip cells presented in 
Figure 1, where the alternate periodic potentials are applied to second, fourth and fifth strips (and 
periodically in other cells), and the excited current of the remaining strips makes the measured signal. 
The system resemble a 3-phase interdigital transducer [6] and conveniently separates the generation and 
the signal circuits. 

Due to the system exact periodicity and different voltages within periodic cells, the summation over 
the repeating strips (over indices k, I) in Equation (14) results in spectral "spikes" at r = mK/M, where 
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M is the number of strips within cells, in the discussed case M — 5 and m — 0, 1, 2, 3, 4. The strip 
voltages and currents can thus be presented in the form (see Figure 1): 

y = a m e~ jimK/A ' mA \ 
h = juje e aS{mK/M)e-^ K l M ^\ (15) 

where /, k = 0,1, 2, 3, 4 and a m are coefficients to be evaluated from the circuit equations: 

V 0 = V 2 = 0, V 1 -V z = U, V 4 = V 3 , I 1 + I 3 + h = 0; 

we intend to measure the combined current I 0 + 1 2 . Naturally, the measured current in real system would 
be the sum over all cells multiplied by the system aperture-width. For convenience of further discussion, 
we define the observed signal as S = (Io + I-zj/u — C, where C is intended to reduce the signal resulting 
simply from the capacitance of the signal strips to the voltage- supplied strips (this reduction can be 
achieved by applying the trimmed capacitance to the supplied strip as shown in Figure 1). 

In the case of finite number of cells, the Equation (14) must be applied with proper summation over 
the applied strip potentials and currents, within cells and over all cells in the sensor. Concerning strip 
potentials, the condition of equal but opposite bus-bar currents feeding the generating strips allows one 
to evaluate the voltage distributions: V ± of upper and lower bus-bars, respectively; V + — V~ — U is the 
applied voltage to bus-bars. This corresponds to superposition of strip voltages = ±U/2 and certain 
bias voltage U°. Noticing that U° would excite almost exclusively the thickness vibration of piezoelectric 
layer and normal longitudinal waves in the body, we may neglect U° yielding only certain signal bias 
like that caused by the strip mutual capacitances. Now, the above-mentioned double summation over 
strip cells results in multiplication of analytic R(r) in Equation (14) by 

L L 

Yl e~ jrin - m)MA = (Lsinx/x) 2 , x = irMLr/K, (16) 

n=l m=l 

where L is the number of strip cells in the sensor. 

The resulting signal S(f) as Figure 3 shows is somehow distorted in comparison with that shown in 
Figure 2(b), but still retains its useful features: max ImlS 1 } appears at frequency close to f a (that is at 
zero frequency deviation in the figure) even for large Im{A; t } yielding good estimation to Re{k t }, and 
the slope of RelS*} at f„ dependence on lm{k t } can be used for estimation of the shear wave damping. 

6. Multiperiodic Electrode System 

In more general case, there are different strips within periodic cells. In the above-discussed sensor, for 
example, strips number 3 and 4 can be joined without any space between them. Such system is no longer 
simple-periodic; it is called [7] "multiperiodic" and the "template electric field" presented above must 
be suitably corrected. Appendix B presents the practical result of rather lengthy analysis [8,9] in form 
of a Matlab® code for direct evaluation of the spatial Fourier expansion coefficients of the corrected 
template function F n , n = 0, 1, 2, 3, . . . (F_ n = F*) for periodic cells of period A. Each cells include 
even number M of strips and they are identically "wired" (by external connections) within all cells. 
This yields identical circuit equations except that the applied voltages can be arbitrarily scaled from 
cell to cell. 
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Figure 3. Example S(f) measured by sensor of 20 cells for several values of lm{kt}; the positions 
of max ImlS 1 } are near zero frequency deviation in all cases. 
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The electric fields is constructed using these Fourier coefficients analogously to Equations (5) 
and (10): 

oo 

n=-oo m (17) 

E(x^ ^ ^ ^ ^ j S n — m oc m F n — m G ^ , 

n=— oo m 

where spectral variable r £ (0, K), K = 2%/ A. There is certain difficulty in numerical evaluation of 
potentials and charges of strips resulting from the above equation. Namely, integration of and 
E(x) may yield infinity if r = 0 (which happens at i — 0). 

This can be avoided well by applying small r instead of 0 when applying the FFT procedure for 
integration (r ~ 1CT 4 is a good choice). Taking the output of FFT{[jS , n _ m F n _ m /(r + nK)]} at 
the given strip center yields the sought potential of the strips in the cell, and taking the difference of 
FFT{ [jF n _ m / (r+nK)] } at the space centers after and before the strip yields the given strip charge, 
and Q\ m ^ respectively, / = 0, . . . , M (expressions in brackets are the input vectors of FFT algorithm). 
These values are exploited in the circuit equations for strip currents and voltages in order to evaluate a m ; 
other necessary equations result from Equation (4) (put p = r + nK) resulting in Equation (11) where 
P n -m should be replaced by F n _ m . 

7. Conclusions 

Modern electronic technology allows one to easily fabricate planar system of strips, which contributes 
to their wide applications in many electronic devices including sensors (SAW gas sensors, for instance) 
and actuators (piezoelectric linear motors, for instance). Rigorous electric field analysis is usually 
necessary for the design and evaluation of electric properties of such devices. In this paper, a method 
of analysis of periodic strips or periodic groups of strips (with arbitrary width and spacing within 
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cells) is presented in interesting application for surface elastographic sensor, which may contribute to 
better accuracy of measurements of tissue properties required in medical investigation and diagnosis of 
skin [10] for instance. 

Acknowledgments 

This work was supported by the Polish Ministry of Science and Higher Education (Grant 
NO. N515 500540). 

References 

1. Erdelyi, E.; Magnus, W.; Oberhettinger, E; Tricomi, EG. Higher Transcendental Functions; 
Erdelyi, E., Ed.; McGraw Hill: New York, NY, USA, 1953; Volume 1, Chapter 3, pp. 166-167. 

2. Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P Numerical Recipes: The Art of 
Scientific Computing; Cambridge University Press: New York, NY, USA, 1992; pp. 252-254. 

3. Bl0tekjaer, K.; Ingebrigtsen, K.; Skeie, H. A method for analyzing waves in structures consisting of 
metal strips on dispersive media. IEEE Trans. Electron. Devices 1973, 20, 1 133-1 138. 

4. Danicki, E.J. Spectral Theory of Interdigital Transducers. Avaliable online: http://www.ippt.gov.pl/ 
~edanickiAianickibook.pdf (accessed on 10 August 2012). 

5. Dupps, W.J.; Netto, M.V.; Herekar, S.; Krueger, R.R. Surface wave elastometry of the cornea in 
porcine and human donor eyes. MSEJ Refract. Surg. 2007, 23, 66-75. 

6. Lewis, M. SAW Filters Employing Interdigitated Interdigital Transducers. In Proceedings of the 
HIE Ultrasonic Symposium Proceedings, San Diego, CA, USA, 27-29 October 1982; pp. 12-11. 

7. Aoki, T; Ingebrigtsen, K.A. Acoustic surface waves in split periodic metal gratings on a 
piezoelectric surface. IEEE Trans. Sonic. Ultrason. 1977, 24, 179-193. 

8. Danicki, E.J.; Langli, B.; Bl0tekjaer, K. Spectral theory of EM wave scattering by periodic strips. 
IEEE Trans. Antenna. Propagat. 1995, 43, 97-104. 

9. Tasinkevych, Y. Electrostatics of Planar System of Conducting Strips. In Electrostatics: Theory 
and Applications; Bertrand, C.L., Ed.; Nova Science Publish Inc.: New York, NY, USA, 2011; 
pp. 189-221. 

10. Timofeev, G.A.; Faustova, E.E.; Fedorova, V.N.; Timofeev, A.B.; Kornopeleva, L.S. Investigation 
and comparative characteristic of factors determining human skin elasticity (in Russian). Biofizika 
2008, 53, 133-138. 

Appendix A 

Consider an isotropic elastic halfspace z < 0 characterized by Lame constants A,/i and mass 
density p, resulting in acoustic wave velocities v t = a/ (2/x + X)/p and v t = \fpjp~, respectively for 
longitudinal and shear waves. The plane acoustic wave-fields include displacements u 1>2 and normal 
traction T 3i , i = 1,3. The harmonic wave-field is considered in the form: exp(jut — jpx — jqiz) where 
qi = sj (ui/vi) 2 — p 2 , i = l,t with sign chosen to satisfy the radiation condition [4] at z — > oo; p is 
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arbitrary. The wave-field is the superposition of longitudinal and shear wave-components (we drop the 
time-dependence) : 







P ~Qt 




u 3 




_ Qi P 





e -jqiz 0 



0 



-jqtz 



T 31 
T 33 



^u 3) i + nu 1)3 
Ami,! + (A + 2fi)u 3>3 





Fi 














f ti = df/dxi 



(18) 



(xi = x,x 3 = z). For sliding contact with the piezoelectric layer at z = 0, one easily obtains the 
p-dependent function evaluated at z = 0: 



Z = T33/M3. 



(19) 



Much more complicated is the characterization [4] of piezoelectric layer; the theory yields the 
dependence of electric tangential field (E = jpip) and normal induction (D = D 3 ) on the layer upper 
side, which is considered traction-free on normal displacement (u 3 ) and traction T 33 applied to its bottom 
side which is considered metalized: 



-3P = y y n • (20) 

(p ^ z 2 i z 2 2 j [ i-> 3 

Substitution of Equation (19) yields Equation (9). 
Appendix B 

Consider even number M of strips placed within the domain [0, 1]: the strip left and right edges 
are specified in two-column matrix g of M rows (for the algorithm convenience the edge positions are 
modified in the line two). The algorithm presented below yields the data vector F = [F n ], for specified 
numbers of harmonics: n = 0 N — 1. 



ne=length (g) ; 

g=g-sum ( sum (g' ) ' /2 ) /ne; 

h=li*pi*sum (g' ) ' ; 

y= [ li; zeros (N-l , 1 ) ] ; 

for m=l : ne ; 

t=cos (pi* (g (m, 2) -g (m, 1) ) ) ; 

a=li; 

b=li*t; 

z= [ zeros ( 1 , N-2 ) b*exp (h (m) 
for 1=2: N-l; 

c= (t* (2*1-1) *b- (1-1) *a) /l; 

a=b; 

b=c ; 

z (N-l) =c*exp (l*h (m) ) ; 
end; 



%ne even: 



%shif t ing 
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for i=N:-l:l; 

y (i) =2*z (1+N-i :N) *y (1 :i) ; 

end; 

end; 

y= [ zeros (ne/2 , 1 ) ; y ] ; 
F=y(l:N) ; 
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